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pL| I 6 Abstract In this paper, the closed-form analytic solutions of two new Faraday's stand- 

' 7 ing solitary waves due to the parametric resonance of liquid in a vessel vibrating ver- 

^ ' 8 tically with a constant frequency are given for the first time. Using a model based on 

9 the symmetry of wave elevation and the linearized Boussinesq equation, we gain the 

closed-form wave elevations of the two kinds of non-monotonically decaying standing 

T3 ; 11 solitary waves with smooth crest and the even or odd symmetry. All of them have 

^ ; 12 never been reported, to the best of our knowledge. Besides, they can well explain some 

C 13 experimental phenomena. All of these are helpful to deepen and enrich our under- 

. 14 standings about standing solitary waves and Faraday's wave. 

15 PACS Number: 47.35.Bb 

^ 16 Key words Standing wave, solitary, Faraday's wave, parametric resonance 
Oh 



> . 

^ IB 1 Introduction 

i> : 

(N . 

' 19 As pointed out by Faraday [T] and Benjamin &; Ursell [2], when a vessel containing 

O 20 liquid vibrates vertically with a constant driving frequency il, the so-called parametric 

21 resonance occurs so that standing surface waves are observed, in case that the liquid 

22 oscillates with a constant frequency u that is half of the driving frequency fi, say, 

23 a; = VL/2. In fact, Faraday waves have been observed in many fields of science. For 

24 example, the experimental observation of Faraday waves in a Bose-Einstein conden- 

25 sate was reported by Engels et al p], and the Faraday instability on a free surface of 

26 superfluid "^He was investigated by Abe et al [1] and Ueda et al |5]. In 2011, using a 

27 vertically vibrating Hele-Shaw cell (i.e. nearly two dimensional) partly filled with wa- 
2B ter, Rajchenbach, Leroux and Clamond [6] did an excellent experiment and observed 

29 two new standing solitary surface waves with the odd or even symmetry. These new 

30 standing waves have an unusual characteristic: their elevations non-monotonically 

31 decay to zero in the horizontal direction, while vibrating periodically in the vertical 

32 direction. Especially, they pointed out that "the existence of an oscillion of odd parity 

33 had never been reported in any media up to now". To the best of our knowledge, 

34 theoretical solutions have never been found for these new standing solitary waves. 

35 In this paper, the closed-form analytic solutions of two new Faraday's solitary 



*To whom correspondence should be addressed, sjhao@sjtu.edu.cn 



2 



36 waveqll due to the parametric resonance of liquid in a vessel vibrating vertically with 

37 a constant frequency are reported. Using a model based on the symmetry of wave 
3B elevation and the linearized Boussinesq equation [7] , we gain the closed- form solution 

39 of two kinds of non-monotonically decaying standing solitary waves with the even or 

40 odd symmetry. Both of them have never been reported, to the best of our knowledge. 

41 Especially, they can well explain some experimental phenomena currently reported by 

42 Rajchenbach, Leroux and Clamond [6]. All of these are helpful to deepen and enrich 

43 our understandings about standing solitary waves and Faraday's wave. 



44 2 Closed-form solutions of the new standing waves 

45 Consider a two-dimensional (2D) Faraday's waves in the water depth /i, excited by 

46 a vertically vibrating horizontal bottom being purely sinusoidal with a single driving 

47 frequency f2. Let u denote the frequency of the excited standing wave, respectively. 
4B In theory, it is well-known that the parametric resonance occurs when u = fi/2, 

49 say, the driving frequency Q of the bottom is twice the frequency u of the liquid 

50 vibration [Il[2], which was currently confirmed once again by the excellent experiment 

51 of Rajchenbach, Leroux and Clamond f6]. Thus, in this paper, we focus on the case 

52 u = Q/2 for the parametric resonance. 

53 We use here a model based on the symmetry of wave elevation and the linearized 

54 Boussinesq equation [7]. Let ri{x,T) denote the dimensionless wave elevation, where 

55 T = u t = Q t/2 denote the dimensionless time and x is the dimensionless horizontal 

56 coordinate with x = corresponding to the wave crest. According to the excellent 

57 experiment of Rajchenbach, Leroux and Clamond [6], we have many reasons to assume 

58 that the wave elevation has either the even symmetry about the wave crest x = 0, i.e. 

ri{x,T) = ri{—x,T), — oo < X < +00, (1) 

59 or the odd symmetry 

ri{x,T) = ~ri{—x,T), — oo < X < +00, (2) 

60 respectively. Assuming that the crest is smooth, the even symmetry ([T]) gives us the 

61 boundary condition 

^.(0,r) = 0. (3) 

62 Besides, the odd symmetry ([2]) is equivalent to the boundary condition 

v{0,t) = 0. (4) 

63 Using the above symmetry and the boundary condition at x = 0, we only need seek 

64 a solution ri{x, r) in the interval < x < +oo. It should be emphasized that the 

65 symmetry plays an important role in our approach, as shown below. 



' Here, the solitary wave means the localized wave. 
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66 Assumed that the fluid is inviscid, incompressible, and the flow is irrotational in 

67 < X < +00 (i.e. the flow is not necessarily irrotational at x = ). Such kind of free 
6B surfaces in the interval < x < +oo can be modeled approximately by the famous 

69 Boussinesq equation [7] , which describes many wave phenomena in shallow water. In 

70 physics, the principle of relativity requires that the equations describing the laws of 

71 physics have the same form in all admissible frames of reference. Therefore, following 

72 Boussinesq [7] and using water depth h as a. characteristic length, one can gain the 

73 dimensionless Boussinesq equation in the reference-frame fixed with the vertically 

74 vibrating bottom: 

Vtt - g' (^Vxx + ^Vxxxx + ^VVxx + ^VxVx^ =0, < X < +00, (5) 

75 subject to the bounded condition 

\ri{x,T)\<S, 0<x<+oo, (6) 

76 where g' is the so-called dimensionless "apparent gravity acceleration" and the is a 

77 large enough positive constant, respectively. Note that the above equation has exactly 
7B the same form as the traditional Boussinesq equation [7], except that the "gravity ac- 

79 celeration" term g' has a different meaning. Obviously, according to Einstein's theory 

80 of general relativity, for an observer moving with the vertically vibrating horizon- 

81 tal bottom that is not an inertial frame of reference, the so-called apparent gravity 

82 acceleration reads 

g' = G{1 - Fcos2r) , 

83 where G = g/{hu'^) is the dimensionless acceleration of gravity, g is the acceleration 

84 due to gravity, F = T/g denotes the dimensionless driving acceleration with F being 

85 the amplitude of the forcing acceleration of the bottom, respectively. 

86 Assume that the wave amplitude is so small that all nonlinear terms of ([5]) can 

87 be neglected. Thus, we have the linearized Boussinesq equation in the non-inertial 

88 frame of reference fixed with the vertically vibrating bottom: 

,..-G(l-Fcos2r)(,.. + l,„)=0, < . < +00. (7) 

89 Our purpose is to find the solutions of Eq. ([7]), subject to the bounded condition (Q 

90 and either the boundary condition ([3]) for the standing solitary waves with the even 

91 symmetry of elevation or (jl]) with the odd symmetry, which oscillate periodically in 

92 time r with the period T = 27r. 

93 Note that t = ut and u = Q/2, where Q is the driving frequency of the verti- 

94 cally vibrating bottom. According to the excellent experiment done by Rajchenbach, 

95 Leroux and Clamond [6], the parametric resonance occurs when u = Q/2. Thus, we 

96 express the standing wave in the form 

ri{x, r) = /(r) e^"", < x < +oo, (8) 
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97 where /(r) is a periodic function with the period T = 2tt, and A is an unknown 

9B eigenvalue to be determined. Note that the eigenvalue A can be real or complex. 

99 Substituting the above expression into ([7]), we have a linear ordinary different equation 

fir) - GX' (l + (1 - Fcos2r)/(r) = 0. (9) 

100 The above equation can be rewritten as the standard Mathieu equation 

r(r) + [a-2gcos(2r)]/(r) = 0, (10) 

101 where 

a = -G\' (l + y) , <1= (0 «• (11) 

102 Thus, /(r) is a periodic Mathieu function with the characteristic value a and the 

103 parameter q, denoted by /(r) = Mc{t; a, q). Note that similar Mathieu-type analyses 

104 have been carried out for Faraday waves in one and two-component Bose-Einstein 

105 condensates P-[TU] 

106 It is well-known that, for a given non-zero parameter g, the corresponding Math- 

107 ieu functions /(r) of Eq. ffTOj) are periodic in r only for certain values of a, called 

108 Mathieu characteristic values. According to Floquet's Theorem, any Mathieu func- 

109 tion /(r) can be written in the form e*'"^/*(r), where /*(t) has period 27r and r is the 
no Mathieu characteristic exponent. The Mathieu function /(r) is periodic only when 

111 the characteristic exponent r is an integer or rational number. 

112 For given characteristic value a and parameter q = {aF)/2, let us consider the 

113 even Mathieu function /(r) of ffTOj) with the characteristic exponent r = 1 so that 

114 /(r) has the period T = 27r. As mentioned above, the characteristic value of the even 

115 Mathieu function /(r) with characteristic exponent r = 1 given by the parameter q 

116 must be equal to the characteristic value a itself. This gives, by means of the computer 

117 algebra system Mathematica, the following nonlinear algebraic equation 

MathieuCharacteristicA[l, aF/2] = a, (12) 

118 where the Mathematica command MathieuCharacteristicA[r, q] is used to gain the 

119 characteristic value a for even Mathieu functions with characteristic exponent r and 

120 the given parameter q. Given the dimensionless driving acceleration F, the above 

121 nonlinear algebraic equation contains only the unknown characteristic value a, denoted 

122 by a*. The corresponding solution 

/(r) = M,(r;a*,Fa72) 

123 is an even Mathieu function with the period 27r. It should be emphasized that the 

124 characteristic value a* depends only on the dimensionless driving acceleration F. 

125 It is found that the nonlinear algebraic equation (fT2l) has two solutions in general. 

126 For example, when F = 2, we have a positive characteristic value a* = 2.49527 and a 

127 negative characteristic value a* = —3.47044, respectively. For different values of the 
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dimensionless driving acceleration F, we have different characteristic value a*. The 
two curves of the characteristic value a* versus the dimensionless driving acceleration 
F are as shown in Figs. [TJandO Note that, the maximum positive characteristic value 
a^^^. is 2.52168, corresponding to the dimensionless driving acceleration F = 2.28. 
Note that, according to the excellent experiment of Rajchenbach, Leroux and Clamond 
|6], the parametric resonance was found when F = 2.0. So, the above theoretical 
result can partly explain why Rajchenbach, Leroux and Clamond |E] observed the 
parametric resonance in case of F = 2.0, since the corresponding characteristic value 
a* = 2.49527 is rather close to a^^^. = 2.52168 so that the parametric resonance is 
easier to be created and observed. We will discuss this later in details. 



Figure 1: The positive characteristic values of a* versus the dimensionless driving 
acceleration F given by the linearized Boussinesq equation ([7]). 

13B According to f lTT]) . as long as the characteristic value a* is known, it is easy to 
139 gain the unknown eigenvalue A by solving the nonlinear algebraic equation 



AM1 + 



+ fi = 0, 



(13) 



where 



141 Thus, A is dependent upon the dimensionless gravity acceleration G and the dimen- 

142 sionless driving acceleration F, since a* is determined by F only. When /i > 3/4, the 

143 above nonlinear algebraic equation have four complex roots 



2,3,4 



± 



-- ± iW3 f u — - 
2 V V 4 



(14) 
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Figure 2: The negative characteristic values of a* versus the dimensionless driving 
acceleration F given by the linearized Boussinesq equation ([7]). 

144 where i = ^/—l denotes the imaginary unit. When < /i < 3/4, there exist four pure 

145 imaginary roots 



Ai,2 = ±i 

146 When /i < 0, there are two pure imaginary roots 



1 




(15) 



A 



3,4 



±i 



\ 




(16) 



and two real roots 



A 



3,4 



± 



\ 




(17) 



14B Thus, when /i > 3/4, we have four complex roots A = ±a ± /3i with a > and 
149 /3 > 0, corresponding to the wave elevation in a general form 

r]{x,T) = Ai Mc{r; a* ,q*)e~°'^ {cos f3x + ism (3 x) 

+ A2 Mc{T;a*,q*)e~'^'' {cos (3x - ism (3x) 

+ A3 Mc{t; a*, g*)e+"'' (cos /3x + i sin /3x) 
+ Mc(r; a*, g*)e+"'^ (cos sin (18) 



150 where Ai,A2, A3 and are constants. However, restricted by the bounded condition 
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151 (Q, only the solution in the form 

r]{x,T) = Ai Mc{t] a* ,q*)e~°'^ {co?, fix + ism (5x) 
+A2 Mc{t; a*, q*)e~°"' (cos f3x - i sin (3x) 
= Mc{t; a*, q*)e~°''' (Aq cos (3x + Bq sin I3x) , < x < +00, (19) 

152 has physical meanings, where Aq = Ai + A2 and Bo = {Ai — A2)i are real constants. 

153 Thus, using the boundary condition (jl]) for the odd symmetry and enforcing 

154 Aq = 0, we have the wave elevation 

ri{x, t) = Bo Mc{t; a*, q*) sin(/3x)e-"'^, < x < +00. (20) 

155 Then, due to the odd symmetry ([2]), we have the odd-pattern elevation 

?7(x, r) = 5o ^c(r; a*, g*)sin(/3x)e""l''l, -00 < x < +cx). (21) 

156 Similarly, using the boundary condition ([3]) for the even symmetry ([T]), we have the 

157 even-pattern wave elevation 

rj{x,T) = Ao Mc{T;a*,q*)e'^ 

158 which has a smooth crest and is valid in the whole domain — 00 < x < +00. Note 

159 that the standing solitary wave elevations 0211) and fl2^ decay non-monotonically in 

160 the X direction, and have no peaked crest. 

161 For example, in case of the driving frequency ri = 11 Hz with the vibration 

162 amplitude 4.1 mm and water depth 5 cm, which were used by Rajchenbach, Leroux 

163 and Clamond [6j in their excellent experiment, we have the dimensionless driving 

164 acceleration F ^ 2 and the dimensionless gravity acceleration G ~ 0.164. When F = 

165 2, there exist one positive characteristic a* = 2.49527 and one negative characteristic 

166 a* = —3.47044, corresponding to fi = a*/G = 15.2151 and fi = —21.1612, respectively. 

167 Especially, when a* = 2.49527, i.e. fi = 15.2151, we have four complex eigenvalues 

A = ±1.62113 ±2.03126i, 

168 corresponding to a non-monotonically decaying standing solitary wave with the odd 

169 symmetry and the smooth crest 

7]{x,t) = AoMe(r; 2.49527, 2.49527) sin(2.03126x) e-^-^2"3|^l, (23) 

170 as shown in Fig. [3l and the wave elevation with even symmetry and smooth crest 

r/(x,r) = AoM,(r; 2.49527, 2.49527) e"^-^^"^l^'l 

X [cos (2.03126x) + 0.798091 sin (2.03126|x|)] , (24) 

171 as shown in Fig. HJ respectively. 



a 



cos{(3x) + ( ^ ) sin(/3|x|) 



(22) 
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Figure 3: Non-monotonically decaying standing solitary wave ( 123|) with the odd sym- 
metry and a smooth crest, (a): r = 0; (b): r = 7/5; (c): r = 9/5; (d): r = tt. 
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Figure 4: Non-monotonically decaying standing solitary wave f l2^ with the even 
symmetry and a smooth crest. Sohd hne: r = 0; Dashed hne: r = 7/5; Dash-dotted 
hne: r = 37r/4; Dash-dot-dotted hne: r = tt. 

172 It should be emphasized that the standing solitary wave (1231) has the odd parity 

173 about a; = 0. Note that Rajchenbach, Leroux and Clamond [6j found a similar 

174 standing solitary wave with the odd parity in their excellent experiment, and pointed 

175 out that "the existence of an oscillon of odd parity had never been reported in any 

176 media up to now". Thus, the closed-form solution (123|) might provide a theoretical 

177 explanation for this experimental phenomenon. 

178 Note that the standing solitary waves fl23|) and fl24|) do not decay monotonically, 

179 as shown in Figs. [3] and HI which are qualitatively similar to those experimentally 

180 found by Rajchenbach, Leroux and Clamond [6J. Note that these non-monotonically 

181 decaying standing solitary waves ( l23l) and ( 12^ are not exactly the same as those 

182 found by the excellent experiment of Rajchenbach, Leroux and Clamond [6]. Such 

183 kind of difference may likely be attributed to the probe motion in their experiment [6] , 

184 and also to the neglect of the nonlinearity of the Boussinesq equation that is valid 

185 for fairly long waves with small-amplitude in shallow water. The nonlinearity of the 

186 Boussinesq equation might affect the eigenvalue A, which determines the decay-rate 

187 of wave elevation. However, the nonlinear terms should not qualitatively influence the 

188 shape of wave elevation. Since we mainly focus on the shape of wave elevation in this 

189 article, the neglect of the nonlinear terms is acceptable. Obviously, better analytic 

190 approximations of the two new standing solitary waves should be gained, if the exact 

191 Boussinesq equation ([5]) or the fully nonlinear wave equation is solved. 
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192 3 Discussions and concluding remarks 

193 In this paper, some new Faraday's waves due to the parametric resonance of hquid in 

194 a vessel vibrating vertically with a constant frequency are reported. Using a model 

195 based on the symmetry of wave elevation and the linearized Boussinesq equation, we 

196 gain the closed-form solutions of two kinds of non-monotonically decaying standing 

197 solitary waves with the even or odd symmetry. They can well explain, although 

198 partly, some experimental phenomena currently reported by Rajchenbach, Leroux 

199 and Clamond [B]. 

200 First, our closed-form solution ( ETj) has the odd parity about x = 0, as shown 

201 in Fig. 131 Note that Rajchenbach, Leroux and Clamond [6] found a similar standing 

202 solitary wave with the odd parity in their excellent experiment, and pointed out that 

203 "the existence of an oscillon of odd parity had never been reported in any media up 

204 to now" . So, our closed-form solution fl2Tl) supports this experimental phenomenon. 

205 Secondly, based on the linearized Boussinesq equation, the characteristic value a* 

206 is dependent upon the dimensionless driving acceleration F only. Thus, for given di- 

207 mensionless gravity acceleration G = g/{huj'^) = Ag/ihVt^), the occurrence of the non- 
20B monotonically decaying standing solitary waves mainly depends on the dimensionless 

209 driving acceleration F = T/g of the vertically vibrating bottom: the larger a*, the 

210 larger possibility of the occurrence of the non-monotonically decaying standing solitary 

211 waves, since fi = a*/G > 3/4 is the criterion for the linearized Boussinesq equation. 

212 According to Fig. (H the non-monotonically decaying standing solitary waves occur 

213 with the maximum possibility at F ^ 2.28, corresponding to the maximum character- 

214 istic value a'^^x — 2.52168. Note that Rajchenbach, Leroux and Clamond ||6j observed 

215 the two non-monotonically decaying standing solitary waves at F = 2, corresponding 

216 to the characteristic value a* = 2.49527 that is rather close to a'^ax ~ 2.52168. So, 

217 our theoretical result can explain this experimental phenomena quite well. 

21B Thirdly, Rajchenbach, Leroux and Clamond [6] found experimentally that the two 

219 non-monotonically decaying standing solitary waves occur in an interval Fl < F < 

220 Fji. They gave it a theoretical explanation using Meron's stability theory pQj. Based 

221 on the linearized Boussinesq equation, we gain the criterion of occurrence of the two 

222 non-monotonically decaying standing solitary waves: fi = a* /G > 3/A. According to 

223 Fig. [H a* has a maximum a'^^^ = 2.52168 at F = 2.28. So, given a proper value of G, 

224 one might find a closed interval of F for the occurrence of the two non-monotonically 

225 decaying standing solitary waves. Thus, our theoretical result about the criterion 

226 yU > 3/4 can explain this experimental phenomenon, too. 

227 Seriously speaking, profile of the standing waves should be dependent upon not 
22B only the dimensionless accelerations G, F but also the wave height. However, based 

229 on the linearized Boussinesq equation, the profile of the two standing waves is depen- 

230 dent on G and F only. This is similar to the periodic travelling waves, whose wave 

231 profile is sinusoidal and independent of wave height. Besides, the detailed evolution 

232 of the standing waves should be also influenced by wave height. So, if the influence 

233 of wave height is considered, the nonlinear Boussinesq equation or even the fully non- 
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234 linear water wave equations should be used for more accurate wave profile and better 

235 understandings about such kind of standing waves. 

236 Note that, using the even or odd symmetry, we have either the boundary condition 

237 ([3]) or (jlj) at X = 0, so that it is enough for the governing equation to be satisfied in 

238 the interval < x < +oo except x = 0. This is well-known and widely used in the 

239 theory of differential equations. 

240 Traditionally, one need give a global expression of a solution in the whole domain. 

241 However, this is difficult in many cases. Fortunately, this traditional idea is out of date. 

242 In modern mathematics, we often express a smooth function by a lots of local simple 

243 functions in a finite number of sub-domains: this idea is widely used in Finite Element 

244 Method (FEM). Although there exists singularity at each boundary of the sub-domain 

245 where the solution is not smooth, such kind of solutions are widely accepted and 

246 used. In this paper, we use different base functions to express the solution of a new 

247 type of standing waves found by experiment in two sub-domains only, i.e. (— oo, 0) 

248 and (0, +oo), and connect them by the symmetry and the smoothness condition at 

249 X = 0. Comparing to the FEM, the mathematical approach used in this article is 

250 more traditional: we use the symmetry to divide the whole domain into only two 

251 sub-domains, and besides the solutions are smooth at x = 0. 

252 Finally, it should be emphasized that the closed-form wave elevations of the two 

253 non-monotonically decaying standing solitary waves (!2T|) and (!22l) are obtained under 

254 the assumption of the even or odd symmetry of wave elevation by means of the 

255 linearized Boussinesq equation with the neglect of viscocity of fluid in the interval 

256 < X < -|-oo. The symmetry has an important role in our approach. The fact that the 

257 two closed-form solutions well explain some phenomenon of the excellent experiment 
25B of Rajchenbach et al [6] indicates the validity of this model. In addition, since the 

259 linearized Boussinesq equation is only a simplified model for shallow water waves, 

260 all conclusions and theoretical predictions reported in this article should be further 

261 checked and verified by fine numerical simulations and physical experiments in future, 

262 even though our closed-form solutions well explain some experimental phenomena of 

263 Rajchenbach et al [6]. All of these are helpful to deepen and enrich our understandings 

264 about standing solitary waves and Faraday's wave. 
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